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Abstract 

In supersymmetric theories, the presence of axions usually implies the existence of a 
non-compact, (pseudo)moduli space. In gauge mediated models, the axion would seem a 
particularly promising dark matter candidate. The cosmology of the moduli then constrains 
the gravitino mass and the axion decay constant; the former can't be much below 10 MeV; 
the latter can't be much larger than 10 13 GeV. Axinos, when identifiable, are typically heavy 
and do not play an important role in cosmology. 



1 Introduction 



The strong CP problem, as presently understood, requires either an axion, a massless u quark, 
or spontaneous breaking of CP at some high energy scale. Lattice gauge calculations appear 
to rule out a massless u quark[l]; at least in the widely explored landscape framework, the 
spontaneous CP solution seems unpromising [2]. So in many ways, the axion solution to the 
strong CP problem seems the most likely of the three to be realized in nature. From astrophysics, 
there is a lower bound on f a of order 10 9 GeV. If cosmology is conventional (in particular, if 
the universe was in thermal equilibrium up to some very high temperature, and not allowing 
for possible anthropic selection effects), the axion decay constant is less than about 10 13 GeV. 
In unconventional cosmologies [3, 4, 5, 6], or allowing for anthropic selection[7, 8, 9], it may be 
larger. 

One can attempt to implement the Peccei-Quinn (PQ) solution of the strong CP problem 
in a variety of frameworks: string theory, low energy effective field theory, and models with or 
without supersymmetry. Some of the issues involved in each of these possibilities are surveyed 
in [2]. In this paper, we focus on models of low energy supersymmetry, with gauge interac- 
tions as the messengers of supersymmetry breaking. Unlike intermediate scale ("super gravity") 
breaking, there is no particular compelling dark matter candidate, so the axion seems worthy 
of study. In this context, the scale of PQ breaking, f a , is necessarily far larger than the scale 
of supersymmetry breaking. So the axion lies in a supermultiplet, with a scalar field known as 
the saxion, and a fermion, the axino. As we will see, there is always a modulus in such models 
which determines the value of f a , but this modulus need not be the saxion. 

In gauge mediated models, theoretical considerations constrain the underlying scale of 
supersymmetry breaking, to lie roug hly between about 10 5 GeV and 10 9 GeV. 1 The 

light modulus associated with PQ symmetry breaking is subject to significant cosmological 
constraints[12, 13]. In this paper we will argue, first, that this particle is unacceptable unless 
the axion multiplet couples directly to the messengers and/or the sector responsible for super- 
symmetry breaking. If the axion is to be the dark matter, even this is unacceptable unless the 
scale of supersymmetry breaking lies to the high end of the allowed range, 10 8 GeV or so. 2 

In the rest of this paper, we will explore models in which susy-breaking dynamics and the 

1 Onc will sometimes hear it said that \f\F\ as large as 10 11 is allowed. We are assuming here that the 
underlying microscopic theory does not even approximately conserve CP, or any approximate flavor symmetries. 

2 Scales similar to those considered here are an integral part of the F Theory GUT program of [10, 11]. The 
underlying pictures, however, are quite different, so this convergence of scales is remarkable. 
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couplings to messengers determine f a and the mass of the superpartners of the axion. We will 
construct a variety of models which illustrate the issues we have outlined above. In one class, 
the saxion is a pseudomodulus; in another, the saxion is relatively heavy, with another field 
determining the value of f a . In these models, we can explore the cosmological issues in a sharp 
fashion. Another issue which will concern us is the origin of the hierarchy between the PQ 
scale and the scale of susy breaking. We will explore whether this "PQ hierarchy" might arise 
through small couplings or through tuning, or dynamically. 

The rest of this paper is organized as follows. In section 2, we note that gauge mediated 
models typically possess an approximate R symmetry, which must be spontaneously broken. 
It is natural to ask whether the Goldstone boson of that symmetry, the i?-axion, might be the 
QCD axion. We explain why the answer is almost certainly no. 

In section 3, we present models in which the Peccei-Quinn symmetry is broken at tree 
level, but there is a pseudomoduli space on which f a varies. In one example, the saxion is itself 
a pseudomodulus. In a second, the "i?-saxion" (the light modulus in the Goldstino multiplet) 
is the pseudomodulus which determines /„. The axino and saxion are already massive at tree 
level. At one loop, the mass of the pseudomodulus is estimated. In section 4, we discuss the 
cosmological constraints on gauge mediated models under the assumption that the axion is the 
dark matter. The most severe arises from coherent production of the lightest pseudomodulus. 
We demonstrate that in a generic inflationary model, the saxion starts a distance of order f a 
from its minimum. As a result, it comes to dominate the energy density of the universe around 
the time of nucleosynthesis, with problematic consequences for cosmology, unless its mass is 
large enough. This translates, in turn, into a lower bound on the supersymmetry breaking scale 
as a function of f a . We argue that, as a result, in such a framework, messenger masses are likely 
to be of order f a . In a broad range of circumstances, the axino is short lived on cosmological 
time scales, and does not pose additional problems. Moreover, from this consideration alone, 
f a cannot be too large; / a 's significantly larger than 10 13 GeV seem unacceptable. If the axion 
is not dark matter, the cosmological constraint arising from pseudomoduli is weaker or even 
non-existent. 

In the models of section 3, the hierarchy between the supersymmetry-breaking scale and 
f a is simply put in by hand, the result of a large ratio of masses in the underlying lagrangian. 
One might hope to avoid this in either of two ways. First, the parameters themselves might 
arise through "retrofitting", i.e. they might be dynamically determined by other interactions. 
Alternatively, and more economically, given that the saxion is a pseudomodulus, with a potential 
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varying only logarithmically over large regions of the field space, it would seem that such a 
potential could naturally have a minimum at a hierarchically large value of f a - We explore 
this possibility in section 5, and find that this can be achieved only in a restricted set of 
circumstances. 

Lurking in the background in all of this discussion is the question of whether unknown high 
energy effects spoil the Peccei-Quinn solution of the strong CP problem. As is well-known, this 
question can be organized in terms of higher dimension, CP-violating operators[14]. Even with 
a relatively low scale of super symmetry-breaking and Peccei-Quinn symmetry breaking, it is 
necessary to suppress operators up to rather high dimensions. Assuming that the Peccei-Quinn 
symmetry is an accidental consequence of a discrete Zn symmetry, for example, would seem to 
require, in the models developed here, a huge value of iV[15, 16, 17, 18]. As discussed in [2], 
imposing a constraint that the axion constitute the dark matter, in these types of models, can 
almost account for the quality of the PQ symmetry. If we do not require the axion be the dark 
matter, the required values of N are more modest, but still rather large and now would simply 
appear a peculiar accident. We discuss this issue in section 6. In our concluding section, we 
summarize what we view as the principle lessons of this work, and ask what features of the 
models considered here are likely to be robust. 

2 The R- axion as the QCD Axion 

A theorem of Nelson and Seiberg[19] asserts that generic theories which break supersymme- 
try possess an approximate R symmetry. In any realistic model, this must be spontaneously 
broken 3 . The corresponding Goldstone boson is conventionally called the .R-axion. By way 
of nomenclature, we will refer to the pseudomodulus that accompanies the i?-axion as the 
"i?-saxion" . 

In any case, it is usually said that the i?-axion gains mass when one couples the system 
to super gravity [20]. This is because it is necessary, in order to (nearly) cancel the cosmological 
constant, to insure that the expectation value of the superpotential takes a certain value of 
order m 3 / 2 Mp. This is typically achieved by including a constant in the superpotential, which 
is then suitably tuned. This constant breaks any would-be (discrete or continuous) R symmetry. 
But one might hope, instead, that the R symmetry is spontaneously broken by some dynamics 4 
3 It is necessary to break the R symmetry to give mass to the gauginos. 

4 In a landscape, there might be a large set of theories with a distribution of (W) accounting for the value of 
the cosmological constant. 
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Upon careful consideration, however, this possibility does not appear promising. The issue 
is a mismatch of scales. In models of dynamical breaking, one might expect that 



(W) = A 3 (F) = A 2 



(1) 



In order to cancel the cosmological constant, one needs, instead a hierarchy - a huge hierarchy 
- between the scale of R breaking and the scale of supersymmetry breaking; the scale of R 
breaking, it would seem, should be of order the Planck scale, so f a ~ 10 18 GeV. It is hard to see 
how this could be achieved consistent with the constraints from cosmology and a high quality 
axion. 

3 Models With Hierarchical Breaking of the Peccei-Quinn Sym- 
metry 

As we build models, we will first impose a global, continuous Peccei-Quinn symmetry. Our 
viewpoint will be that this must eventually be accounted for as an accidental consequence 
of an underlying discrete symmetry (or perhaps a continuous gauge symmetry). One of us 
has discussed this issue in the past [16], and we will comment further on this question in the 
concluding section. 

To outline the basic issues, consider, first, a theory with fields S± carrying PQ charge ±1, 
and a neutral field, %■ Take for the superpotential: 



Supersymmetry is unbroken in this theory. The Peccei-Quinn symmetry is broken, and there 
is a moduli space, which can be understood in terms of the "complexification of the symmetry 
group" : 



W = x(S + S^-^ 2 ). 



(2) 



S± = fie 



,±A/u 



(3) 



f a in this model is a function of A 




(4) 



A is naturally referred to as the axion supermultiplet. 
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3.1 Model In Which the Saxion is a Pseudomodulus 

We want to complicate this model so that supersymmetry is broken, at a scale much less than 
A* ~ fa- We will consider, first, O'Raifeartaigh models, postponing (limited) discussion of 
dynamical breaking of supersymmetry and/or the Peccei-Quinn symmetry until later. There 
are various strategies we might adopt. In our first model, we simply add modules to that of 
eqn. 2, which break supersymmetry, fix A (fa), and act as messengers. We will suppose that, 
in addition to the Peccei-Quinn symmetry, the model possesses an R symmetry, as required by 
the theorem of Nelson and Seiberg[19] to break supersymmetry in a generic fashion. Examples 
of sectors which break supersymmetry and the R symmetry are provided by the models of 
[21, 22, 23]. For definiteness, we take this sector to consist of fields A, with R charge 2, and 
01, (f)-i, 03, where the subscripts denote the R charges of the fields, and with superpotential 

W x = AiX(<M_i -F) + mi<f>l + m 2 <p-i<t> 3 - (5) 

All of the fields in eqn. 5 are taken neutral under the PQ symmetry. As shown in [21], in this 
model, there is a range of parameters for which the R symmetry is broken, and 

(X) = x + 9 2 F. (6) 

with x ~ mi. Note, in particular, that a hierarchy between x and F is possible. 

We will take S± to have R charge 0. In order to fix A in eqn. 3, rather than cou- 
ple S± directly to A, we couple these to another set of fields which do couple to A. The 
resulting model has a two dimensional (pseudo) moduli space, one direction associated with 
breaking the R symmetry, the other with breaking the PQ symmetry. The additional fields 
ai, ai, a,2, a 2 , bi, bi, b 2 , 62 j do not transform under the gauge symmetries. We take the additional 
contribution to the superpotential to be 

W 2 = x(S+S- - fi 2 ) + hS + aia 2 + yS- b\ b 2 + A(ojaj + (7) 

We have not labeled the various independent couplings of A to the a^'s, and 6j's to avoid 
cluttering formulas unnecessarily. Also, to avoid cluttering the formulas, we have not used the 
subscripts to label the R or PQ charges, but this lagrangian is consistent with both symmetries. 
For example, denoting R and PQ charges as (R, PQ) we can take: 

ai(l,-l),ai(-l,l); a 2 (-l, 0), a 2 (l, 0); h(l, 1), -1); b 2 (-l, 0), 5 2 (1, 0). (8) 

The superpotential of eqn. 7 is not the most general consistent with the symmetries. For 
example, we can add a term eXS+S-; this is easily seen not to qualitatively alter the behavior 
of the model discussed below, at least for small e. Similarly for other possible terms. 
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In the model as it stands, the scales f a ~ fi and x ~ m; are independent. When we study 
cosmological issues in section 4, however, we will see that cosmological considerations require 
that they be comparable. We can write for S+ and S-: 



then in the limit of small F and for y/hy^x 1 not too large, the minimum of the one loop 
potential for the saxion (real part of cf)) is at 4> = 0. The saxion acquires a mass 



The axino mass is parametrically lighter, by the square root of a loop factor. We will see 
shortly that this is not generic; the axino can easily be - and arguably typically will be - heavier 
than the lightest pseudomodulus. 

We need one final module; it is important that the Peccei-Quinn symmetry be anomalous, 
and that supersymmetry breaking be transmitted to the fields of the MSSM. We accomplish 
this by coupling X and S to a set of fields filling out two 5 and 5's of the Standard Model: 



The fields qi,q~i,£i,£i act as messengers. The fields q2, etc., gain mass as a result of the S± 
expectation values, giving rise to the coupling of the axion supermultiplet, A, to the gauge 
fields of the Standard Model (MSSM). Within this model, the masses of messengers and the 
scale f a are in principle independent. One can contemplate more intricate messenger sectors, 
implementing general gauge mediation[24, 25]. But the most important point about this struc- 
ture is that, given our argument in the next section that f a ~ x, messenger masses are likely 
to be of order f a , in implementations of the Peccei-Quinn symmetry within gauge mediation. 

To summarize the spectrum, we have two pseudomoduli multiplets in this model. At the 
minimum of the Coleman- Weinberg potential: 

1. Superpartners of ordinary fields of the MSSM have masses given by standard model loop 
factors times F/x, as in conventional gauge mediation. 

2. For f a ~ x the saxion mass is of order 




(9) 




1 F 2 



(10) 



W 3 = + hh) + S+q 2 q 2 + S.t 2 l 2 . 



(11) 




(12) 



i.e. in terms of loop counting, larger than that of the MSSM particles. 
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3. The axino mass is the same loop order as those of the MSSM particles: 

m '^^Ta (13) 

though it is further suppressed if R symmetry breaking is small. If \x\ <C f a , the suppres- 
sion is 0(x/f a ). In this situation, in order to obtain suitable gaugino masses, there must 
be additional messengers which couple to X but not to S. 

4. The "r-axion" is massless, as is the Goldstino. The "r-saxion" gains mass at one loop 
similar to that of the saxion. 



3.2 Models in Which the "R-Saxion" Determines f a 



In the model of the previous section, there were two pseudomoduli, one responsible for deter- 
mining f a , one for determining f r , the decay constant of the "r-saxion". This arose because 
of our "modular" structure: by setting certain couplings to zero, it was possible to decouple 
the sector which broke supersymmetry from the sector which broke the PQ symmetry. In the 
absence of supersymmetry breaking, the broken PQ symmetry implies the existence of a moduli 
space (associated with the "complexification" of the symmetry group); the requirement that 
supersymmetry is broken in the other, decoupled, sector, implies the existence of an r-axion, 
according to the theorem of Nelson and Seiberg. In this section, we construct models in which 
the symmetry-breakings cannot be decoupled. The "r-saxion" in these models determines f a ; 
the saxion is not a pseudomodulus at all, and is parametrically much more massive than the 
r-saxion. The axino is massive already at tree level, so much more massive than the lightest 
pseudomodulus. 

The model, again, contains fields S±,x, with quantum numbers as before. In addition, 
there are fields X, Y with i?-charge 2 and PQ charge ±2, respectively. There are additional 
fields, cf)i,(f)-i with the following (R,PQ) charges: 

^(1,0), 0_!(-l,-2); <M1,0), 0_i(-l,2). (14) 

For the superpotential of the model we take: 

W = x(S+S- - fi 2 ) + - XiS 2 _) + Y(0i0_i - A 2 S+ ) + m<j)l + m$. (15) 

Again, we have set some couplings to one to avoid cluttering formulas; our discussion is readily 
modified if these are allowed to vary. This superpotential is the most general (renormalizable) 
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one consistent with the global symmetries and a discrete Z2 symmetry under which the tilde 
fields change sign. Consider the limit of very small Ai = A2 = A, \x ~ m. In this limit (assuming 
that X, Y ~ m, rh), one linear combination of S± combines with x to form a massive field. The 
other linear combination is fixed, but lighter. Classically, there is a one parameter moduli 
space, satisfying X — Y = 0. The modulus is the field 775 (-^ + Y)] this is the R-saxion. The 
R-axion is the imaginary part of this field, and the Goldstino is the corresponding fermion. For 
large X + Y, the field X — Y is the saxion supermultiplet. Because of the constraint, it is not 
a pseudomodulus, and indeed, the real scalar and the axino are massive already at tree level, 
with mass of order The one loop Coleman- Weinberg calculation yields a stationary point 
for the modulus, X + Y. Its mass is of order 

1 F 2 

m?= 16^ (16) 
where x ~ m, rh, and F ~ Xfi 2 . In this region of the parameter space, the hierarchy between F 
and f a is determined by the small couplings A. 

It is interesting to write the effective field theory for small A, integrating out the massive 
fields. Writing S+ = ix + 8S+, S_ = \i + 5S- the massive fields (mass ~ 11) are 

X, H=-L(6S + + 6S-). (17) 
This leaves the light fields (m = or m ~ A/Lt): 

Z is the Goldstino supermultiplet. Its fermion is the goldstino; the imaginary part of the scalar 
is the i?-axion, and the real part the i2-saxion (the lightest pseudomodulus). To obtain the 
effective lagrangian for the light fields, one needs to solve the heavy field equations of motion. 
Solving the equation dW/dH = yields x = v^AZ + 0(^ 1 ), which substituted back in W 
gives 

W = -X(V2n 2 Z + V2ZA 2 - 2/iAB) + Oiii- 1 ). (19) 

This is an O'Raifeartaigh model. It possesses an R symmetry, but the PQ symmetry is explicitly 
broken in the effective theory. The masslessness of the axion, at the level of the low energy 
lagrangian, appears as a consequence of tuning of the parameters. For small Z, the axion is 
principally A, while for large Z, as expected, it is principally B. There is no light axino. 

For small Z, the saxion is a massive field, with mass of order \ 2 ix 2 . For large Z, the field 
A is heavy, while B is light, with mass of order \ 2 n 4 /Z 2 . The real part of Z, the .R-saxion, is 
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still lighter by a loop factor. B is naturally described as the saxion; it is not a pseudomodulus. 
For intermediate values of Z, the saxion and axion are different linear combinations of B and 
A, and so are not precisely aligned; there is not a sharp notion of what one means the by the 
saxion in these cases. 

To generate the coupling of the axion to FF, it is necessary to couple the fields X, Y to 
messengers. A coupling of the form Xqq, for example, generates a coupling of both the axion 
multiplet and the r-axion multiplet to W%. This is relevant not only for the axion, but for the 
light modulus (moduli). We will discuss the cosmological issues raised by these models in the 
next section. As in conventional gauge mediation, there are one loop diagrams contributing 
mass for gauginos and two loop diagrams contributing mass to squarks and sleptons. 

So far, we have not contemplated coupling X to fields other than messengers, so one might 
think that the quarks and leptons should be neutral under the PQ symmetry. However, once 
one attempts to solve the fi problem, some coupling of messenger fields to fields which couple 
to ordinary fields is required. So generically, one expects that the quarks, leptons and Higgs 
will carry PQ charges. Knowledge of these charges is required to fully determine the couplings 
of the axion to ordinary fields; this bears on the question of axion detection. 



4 Saxion/ Axino/ Axion Cosmology 



A priori, the mass scale of the messengers responsible for transmitting supersymmetry breaking, 
and f a , are independent. We have asserted, however, that in practice the relation between these 
scales is subject to cosmological constraints. In this section, we treat the scale of the messengers 
(and thus the scale of supersymmetry breaking) and f a as independent, and make a preliminary 
examination of the consequences. We will see that over much of the interesting parameter range, 
unless the messenger scale is within two orders of magnitude of f a , the pseudomodulus is very 
long lived and dominates the energy density shortly after nucleosynthesis. 

In models of the type we are considering here, the most serious cosmological issues are 
posed by the lightest modulus, which we will denote by P. Recall that the mass of P is of order 

p / f \ /in 12 cvv\ 

= ^^ T « ~ 10 " GeV (iO^Gev) [—fT- ) <20> 
(mp is the mass of P, as opposed to the Planck mass, which we denote M p ). 

During inflation, there is no reason for the lightest pseudomodulus to sit at the minimum 
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of its potential; indeed, we expect it to sit a distance of order f a from its flat-space, zero 
temperature minimum. For example, in the first of the models we considered, during inflation, 
one expects that the potential contains these additional terms (and many others): 

SV = (a\S + \ 2 + b\S-\ 2 )H 2 (21) 

where H is the Hubble constant during inflation; if H < f a the the potential has a minimum at 

A = P^jln(b/a). (22) 

Similar remarks hold for the second model (in which the saxion is not the pseudomodulus) . In 
the class of models where the saxion is a pseudomodulus, P decays principally to an axino and 
a gravitino; in the class where the lightest modulus is the partner of the Goldstino, P decays 
to a gravitino pair. In either case, the decay width is of order 

r 1 m ? in -34 A0 12 GeV \ 2 / m P \ 3 

After inflation, P will be frozen at this point until H ~ mp, at which time P starts to oscillate 
about its minimum. P's energy density at that time is of order mpf 2 . It constitutes a fraction 
of the energy density of order f 2 /M 2 . It behaves like matter, so its fraction of the energy 
density grows with the scale factor. So, e.g., if f a = 10 12 GeV, it represents a part in 1CP 12 of 
the energy density initially. If its mass is, say, 10 -3 GeV, it starts to oscillate when T ~ 10 7,5 
GeV, and dominates the energy density for T ~ 30 KeV. Its lifetime is of order 10 10 seconds, so 
there is a long period of matter domination before the normal time of matter-radiation equality. 
This is followed by a long period in which the universe is dominated by relativistic gravitinos. 
If the mass is closer to 1 GeV, on the other hand, the decay occurs before nucleosynthesis, just 
when the modulus is coming to dominate the energy density. A significant fraction of the decay 
products will include hadrons, so nucleosynthesis will be problematic. 

So quite generally, we seem to want the saxion mass (for f a ~ 10 12 GeV) to be significantly 
greater than 1 GeV. As we have remarked above, if all of the messengers have masses of order 
f a , the saxion will readily satisfy this bound. The situation is potentially problematic, however, 
if some messengers are light compared to f a . For f a = 10 12 GeV and the loop factor 10 -3 , 
for example, one requires y/F > 10 7 ' 5 GeV. This is not much below the scale required if the 
messenger mass is of order f a {\fF = 10 8 - 5 GeV), so it seems most reasonable to suppose that 
the messengers have scale f a . 

For f a ~ 10 9 GeV, there is no significant cosmological constraint arising from the pseudo- 
modulus. 
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For many of the models we consider here, the axino is quite massive, and not problem- 
atic (or even interesting) cosmologically. This seems likely to be generic if f a is comparable 
to the messenger scale. Conceivably there are situations where the axino mass might be sup- 
pressed, and then the axino may significantly constrain, e.g. the reheating temperature after 
inflation [26]. Under such circumstances, it could constitute a significant fraction of the dark 
matter. 

The cosmology of the pseudomoduli requires a more thorough analysis than we have pre- 
sented here. There may well be windows of low messenger scale which are cosmologically 
allowed. But it would seem from this discussion that the most elegant possibility is that the 
PQ scale is comparable to the messenger mass. This does preclude, however, very low scale 
gauge mediation. 

5 Obtaining a Large Hierarchy Dynamically 

In the model of the previous sections, a large ratio of scales in the underlying lagrangian lead 
to a large hierarchy between the f a and y/~F. This could arise dynamically in a model with 
"retrofitted" parameters [27]. One of us will explore this question elsewhere [28], but one might 
hope that it could arise in a somewhat different way. Given that the saxion is necessarily 
a pseudomodulus, and that its potential varies only logarithmically for large values of the 
modulus, PQ breaking would seem highly susceptible to such large hierarchies. In this section, 
we explore strategies to models of this type. 

5.1 Models without New Gauge Fields 

To build simple models, one might introduce a field, S, carrying PQ charge and coupled to 
messengers. Classically, the potential in the S direction vanishes. The coupling to messengers 
leads to a potential for S at one loop. (Similar models, but in which the S field is not a pseu- 
domodulus, have been considered in proposals to understand the fi term [29, 30]). Specifically, 
we write 

X(X 1 M 1 M 1 + A 2 M 2 M 2 ) + ySM x M 2 (24) 

where X, as before, is a field with (X) = x + 6 2 F, and the M's are messengers, taken to 
transform, say, as 5 + 5 of SU(5). This model admits a Peccei-Quinn symmetry with SU(3) 
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and SU(2) anomalies (under which X transforms as well), so the dynamics responsible for 
supersymmetry breaking must respect this symmetry. One can take as transformation laws: 



S -»• e ia S; M x -»• e~ ia M i; X -> e 4a X, 



(25) 



all other fields being neutral. 

Consider, now, the potential for S. At small S, as discussed in [29, 30], the quadratic terms 
vanish. But we are interested in the behavior at large S. This behavior is logarithmic. If the 
coefficient of the log were negative, and if there were higher dimension terms which respected 
the symmetry, then one might get a hierarchically large breaking of the PQ symmetry. The 
coefficient of the logarithmic term is easily found to be: 



So we do not find the desired behavior. This is readily seen to be quite general, in a theory 
with chiral fields only. The relevant contribution to the potential at second order in F 2 is: 



the computation simplifies, as we are only looking for the logarithmic behavior for large values 
of the fields. One can always define a single field, X, to have an F, component, with other fields, 
<f>i, having zero F components. Then it is straightforward to show that the term in Tr M^M 
proportional to X'X is always positive. Allowing for large x, corresponding to a large breaking 
of the R symmetry, does not alter the situation. Such models do not produce the desired 
hierarchy. The same remains true in the models discussed by Shih[21], in which there are fields 
with R ^ 0, 2, which yield R symmetry breaking. The potential still grows logarithmically for 
large values of the fields. 

As an alternative, one can consider models (again like those of [21] which exhibit runaway 
behavior in certain regions of the field space. One might then hope to stabilize the field 
through higher dimension operators. The difficulty with this approach is that typically the 
higher dimension terms give rise to a supersymmetric minimum. This is not surprising. In 
the runaway directions, supersymmetry becomes better and better as the fields become larger. 
As follows from the analysis of [25], the effective theory with the higher dimension operator 
typically has a superpotential of the type: 



V(S) = 




log(S*S)F*F. 



(26) 




(27) 



1 n + 3 



(28) 



(j) n M™ ' 
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This has a supersymmetric minimum. Stabilization can be achieved in a local minimum if, say, 
the leading allowed operator has the form 

_L_ X<f> n+2 

r + ~W~ [ ' 

for some other field X. 



5.2 Models with Additional Gauge Interactions 

With additional gauge interactions (beyond those of the Standard Model), one might hope to 
find runaway behavior in the Coleman Weinberg calculation. It is well-known that the gauge 
corrections to the potential tend to be negative for large values of charged fields. If some of 
the charged fields also carry Peccei-Quinn charge, one might obtain the desired structure. The 
potential, at large fields, might be stabilized by higher order corrections (understood through 
the renormalization group, as in [31, 32]) or through higher dimension operators in the super- 
potential. We have found that it is possible to realize this possibility, but that the resulting 
models are rather complicated. 

Take a simple example, the model of ref. [22, 23], with fields <p±, Z±, etc. Add fields S± , 
i = 1,2, which are assumed to transform, as well, under a PQ symmetry: 

S f e ±ia sf (30) 

For the superpotential, we take 

W = AZ o (0+</>_ - /i 2 ) + MiZ+<f>- + M 2 Z_0+ (31) 

The S fields do not appear in the superpotential at the renormalizable level. For sufficiently 
large [i, (f>± obtain expectation values. If A is not too small, the fields Zq,Z± obtain vanishing 
expectation values. The fields S± can obtain large expectation values, constrained only by the 
U{1) D term. Adding higher dimension operators of the type 

SW = -^{sfs^T (32) 

respects both the gauge and PQ symmetries, and can stabilize the potential. Both symmetries 
are spontaneously broken by large amounts (proportional to a fractional power of M). There are 
a number of fields, in addition to those in the axion multiplet, with masses of order F/f a , times 
loop factors. It is necessary, of course, to couple the fields to messengers, and to impose, 
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say, discrete symmetries which account for this structure. The models are quite baroque, but 
at least provide an existence proof that this sort of model building is possible. 

An alternative, closer in spirit to Witten's inverted hierarchy [31], invokes a new, non- 
abelian (SU(2)) gauge group, Consider the following superpotential 

W = \X{Tr[<j) + <j)-] - /i 2 ) - mTr[Zcb-] - aYTr[^\] (33) 

The fields cj)±, Z are in the adjoint of SU (2) while A and Y are singlets. The PQ charges of <j>± 
are ±1, while Z and Y have PQ charge 1 and —2 respectively. \i and m are mass parameters. 
X, Y, Z have R— charge 2 while <p± have R— charge. 

At the minimum of the potential all the adjoints fields commute and their vev's can be 
simultaneously diagonalized. SU(2) is broken to £7(1) and the pseudomoduli space is three 
dimensional and includes the axion and the R-axion. Exploiting the R and PQ symmetries to 
set A, Y G R. We have: 

Z = Ya 3 , Y = X^ (34) 
4>- = s<t>+, (f>+ = 7^03 (35) 
where s is determined in function of the parameters of the model as follows: 



9^,2 \2,,2 

s(2a 2 + s 2 X 2 ) = 2a A / (36) 



The potential at the minimum is 



2 _ £ W + 3 f 2 A 2 ) 4 
V -\ F \ ~ 4^ m ( 37 ) 

For large A, the tree level spectrum is comprised of: 

1 . Fields which are mostly combinations of , 0_ with masses of order A A 

2. The rest of the fields charged under the unbroken £7(1) with masses of order gY 

3. Massless particles: three massless scalars: the axion, R-axion and R-saxion; the goldstino; 
one massless gauge boson for the unbroken £7(1) with its gaugino. 

4 

4. All the remaining fields neutral under the unbroken U(l) get masses of order 
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The R-saxion and the U(l) gaugino receive masses by quantum effects which we estimate in 
the following. 

For X 3> /x, m the R-saxion potential dependence on X is obtained by replacing the 
parameters in the tree level potential with their values at the scale X. The RG equations are 
such that for a 2 > |<7 2 the effective potential is increasing. There is an IR fixed point at 
a 2 = ^A 2 = ^g 2 while in the UV a and A are either asymptotically free or both diverge. 
In the latter case by choosing the initial conditions with a 2 < |g 2 and decreasing potential a 
minimum will be generated for larger values of X. 

It is not hard to achieve the required hierarchy between X and \/\F\ with all the couplings 
remaining well in the perturbative regime for the range of X of interest. One example of the 
resulting potential and RG running is shown below for g = 1. The R-saxion typically acquires 
a mass of order jj-^y as does the gaugino for the unbroken U(l). 




Figure 1: Plot for the potential in function of log(X). AV = V(X) — V(X m i n ). Also plotted 

2 \ 2 

are u = v = -r and — as a function of los(X) 

We can now couple the singlet Y to messengers in 5, 5 of SU(5). Provided that the 
coupling to the messengers is not too large, the potential will still develop a minimum. There is 
a small kinetic mixing e arising from 3 loop effects between the hidden unbroken U(l) and the 
EM U(l). As a result all of the fields charged under the former acquire EM charges of order 
e. As all fields charged under the unbroken U(l) have masses of the order of the PQ breaking 
scale f a in order to avoid cosmological constraints [33, 34] we need f a > e 2 10 15 Gev. 

This model indicates that it is possible, in principle, to account for the large value of f a 
relative to y/\F\ through an inverted hierarchy. But the model is complex, and introduces new, 
possibly undesirable, light degrees of freedom (such as the massless gauge bosons of the extra 
U(l)). In the end, models such as those we have discussed which account for the PQ hierarchy 
through small parameters are perhaps more plausible. 
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6 Axion Quality 



The most troubling feature of the axion solution to the strong CP problem is the requirement 
that the PQ symmetry be of extraordinarily high quality. In [2], the quality, Q, was defined by 
writing the axion potential as 



In order that the axion solve the strong CP problem, one requires Q < 10 6 , if f a = 10 GeV. 

Adopting the language of the "landscape", states (or theories) with such a high quality 
PQ symmetry are not likely generic; one can ask what might select for them. Since the laws 
of nuclear physics are hardly sensitive to 9, the only plausible answer we can see is that there 
might be classes of theories (states) for which axions are a generic form of dark matter, and that 
something close to the presently observed dark matter density is a requirement for a hospitable 
universe [7, 8, 9]. In string theory (and/or in higher dimension theories) we understand how a 
non- linearly realized PQ symmetry can arise in the four dimensional theory as an accidental 
consequence of features of the microscopic theory. As explained in reference [2], selecting for 
dark matter in such settings may well account for the requisite Q. 

In gauge mediation the situation is different. As we have stressed here, and as discussed also 
in [2], if gauge mediation is the origin of supersymmetry breaking, and an axion is responsible for 
the absence of strong CP violation, then a sensible cosmology requires that any "stringy moduli" 
be fixed at very high energies (supersymmetrically) , and the Peccei-Quinn symmetry must be 
broken within the low energy field theory, as in the models discussed in this paper. Within field 
theory models, the sensitivity to unknown high energy effects is particularly troubling[14]: in 
some range of energy, the theory must exhibit a linearly realized global symmetry of extremely 
good quality. 

To appreciate the severity of the problem, note that PQ violating operators can arise in 
the superpotential 



V a = Qf*cos(--8 ). 



(38) 



J a 




(39) 



yielding contributions to the potential of order 




(40) 
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In the present paper, we have argued that F ~ 10 17 GeV 2 is a natural scale; correspondingly 
we require n > 13! Ignoring the cosmological issues we have raised, even with f a = 10 9 GeV 
and the smallest F's we can contemplate in gauge mediation, we require n > 9. If the Peccei 
-Quinn symmetry results from, say, an underlying discrete symmetry, we need TV > 9 or 
N > 13 to have a viable axion solution to the strong CP problem. 

In a landscape framework, one might expect that such large discrete symmetries are ex- 
tremely rare[39]. On the other hand, in gauge mediation, the axion is a particularly plausible 
dark matter candidate. Again, as described in [2], selecting for such states gives a very flat 
potential, but not quite flat enough to solve the strong CP problem, except for large /„. For 
f a = 10 12 GeV, for example, selecting for axion dark matter requires n > 11. In other words, 
the dark matter requirement comes close, but whether this is close enough is not clear. (With 
supersymmetry broken at an intermediate scale as in "gravity mediation", dark matter does 
not seem a persuasive criterion, as there are other dark matter candidates which seem less 
expensive. 5 ) 

One might argue that a non-generic superpotential, perhaps plausible due to the non- 
renormalization theorems and to experience with string theory [36], could account for the ab- 
sence of such superpotential couplings. However, significant PQ violating terms are likely to 
appear in the Kahler Potential, 



In our standard case where F ~ 10 17 GeV 2 and f a = 10 12 GeV, we must have n > 8. Looking 
at the Kahler potential we see that this corresponds to a Zn of 8, a modest improvement over 
the superpotential constraint. 

7 Conclusions 

One of the great successes of (critical) string theories is that they yield Peccei-Quinn symme- 
tries which hold to a high degree of accuracy [37, 38]; these PQ symmetries can often be thought 

5 Some discussion of axions in gravity mediation, from a different point of view, appears in [35]. 




(41) 



and hence yield a constraint 




(42) 
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of as an accidental consequence of higher dimensional gauge symmetries. They are good sym- 
metries to all orders of perturbation theory, but fail non-perturbatively. Provided that the high 
energy non-perturbative effects are sufficiently small, these theories seem a suitable setting to 
implement the axion solution of the strong CP problem. As discussed in [2], one faces at least 
two issues with such models. First, f a is likely to be large; one needs to reconcile this with 
cosmological constraints, perhaps along the lines of [3, 4, 5, 6] or [7, 8, 9]. Second, one has to 
ask why, once moduli are fixed, there is an axion sufficiently light (a PQ symmetry of sufficient 
quality) to solve the strong CP problem. As discussed in reference [2], one possibility - perhaps 
the only one - is that there is a selection for axion dark matter. For f a > 10 14 GeV, account 
for the quality of the QCD axion. If low energy supersymmetry plays no role in nature, or with 
intermediate scale breaking supersymmetry breaking, 6 this might provide an adequate under- 
standing of the strong CP problem. In gauge mediation, we have seen that the requirement of 
axionic dark matter selects for a very good PQ symmetry, but perhaps not quite good enough 
(unless, again, f a is uncomfortably large). 

We can ask what features of the model of eqn. 3 might be expected to be generic. Among 
these: 

1. f a not much larger than 10 13 GeV. This seems forced by the challenges of saxion cosmol- 
ogy; otherwise, the decay of the saxion is very late, for any plausible scale of supersym- 
metry breaking. 

2. Messenger masses of order the Peccei-Quinn scale: again, this seems forced by the chal- 
lenges of saxion cosmology. If another solution is found to this problem, these conditions 
might be relaxed. 

3. Axinos with mass typically as large or larger than that of the lightest pseudoscalar: this 
seems generic, as a consequence of the first point, which forces the R symmetry breaking 
scale to be comparable to f a . 

In the end, the axion solution to the strong CP problem, within the framework of gauge 
mediation, seems highly constrained. The Peccei-Quinn and messenger scales are likely to be 
similar. If supersymmetry is discovered at the LHC, and if evidence accumulates for a gauge- 
mediated super symmetric spectrum, this could well point to a detectable axion as the dark 
matter, with the NLSP decaying well outside of the detector. 

6 As we have remarked, in the intermediate scale case, it may be hard to understand why the axion is a more 
generic form of dark matter than the conventional neutralino. 
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